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O ■ ABSTRACT 
O ■ 

<N : 

Q . Recently, the entropy excess detected in the central cores of groups and clus- 

O ! ters has been successfully interpreted as being due to radiative cooling of the 

^ \ hot intragroup/intracluster gas. In such a scenario, the entropy floors Sfloor in 

^ \ groups/clusters at any given redshift are completely determined by the conser- 

^ \ vation of energy. In combination with the equation of hydrostatic equilibrium 

,—1 \ and the universal density profile for dark matter, this allows us to derive the 

^ \ remaining gas distribution of groups and clusters after the cooled material is re- 

T— I \ moved. Together with the Press-Schechter mass function we are able to evaluate 

^ \ effectively how radiative cooling can modify the predictions of SZ cluster counts 

^ \ and power spectrum. It appears that our analytic results are in good agreement 

Ph| with those found by hydrodynamical simulations. Namely, cooling leads to a 

O \ moderate decrease of the predicted SZ cluster counts and power spectrum as 

\ compared with standard scenario. However, without taking into account energy 

\ feedback from star formation which may greatly suppress cooling efficiency, it is 

\ still premature to claim that this modification is significant for the cosmological 

^ \ applications of cluster SZ effect. 
\ 

Subject headings: cosmology: theory — cosmic microwave background - galaxies: 
clusters: general — X-rays: galaxies 



1. Introduction 

With the rapid development of observing techniques (Birkinshaw, Hughes & Arnaud 
1991; Jones et al. 1993; Birkinshaw & Hughes 1994; Carlstrom, Joy & Grego 1996; Myers et 
al. 1997; Hughes & Birkinshaw 1998; etc.), the Sunyaev-Zel'dovich (SZ) effect has become 
one of the most powerful tools for the detections of high redshift clusters (Joy et al. 2001) 
and cosmic microwave background (CMB) anisotropy on small scales (Mason et al. 2002). 
Indeed, the redshift-independence is the major advantage of non-targeted SZ surveys over 
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traditional optical and X-ray observations. This arises from the fact that the SZ effect 
depends uniquely on the intrinsic properties of the warm-hot gas associated with cosmic 
structures, whilst the photons interacting with the gas come from CMB at very high redshift 
z ~ 1000. Because robust constraints on cosmological models are provided by the most 
massive and distant clusters, growing interest over the past years has been focused upon 
how well the fundamental cosmological parameters can be constrained with non-targeted SZ 
cluster surveys (Molnar, Birkinshaw & Mushotzky 2002) and SZ power spectrum (Bond et 
al. 2002). The sensitivity of the expected SZ cluster counts and SZ power spectrum to the 
underlying cosmological model is quite impressive. 

On arcminute scales, the strength of thermal SZ signals is directly proportional to the 
total thermal energy of the hot gas confined in clusters. Unlike the dark matter component 
of clusters whose dynamical behavior is governed purely by gravity, the intracluster gas is 
easily disturbed by many complicated physical processes in addition to gravity and thermal 
pressure, which include (non)gravitational heating, radiative cooling, star formation and 
energy feedback, magnetic fields, etc. The reliability of cosmological applications of SZ 
surveys on arcminute scales is closely connected with the question of how well one can handle 
these non-gravitational mechanisms. At present, a sophisticated treatment of the problem 
must rely on hydrodynamical simulations coupled with some semi-analytic approximations 
(e.g. da Silva et al. 2000, 2001; Seljak, Burwell & Pen 2001; Zhang, Pen & Wang 2002; 
White et al. 2002). Yet, the effective probe of cosmological models using non-targeted SZ 
surveys can be more easily achieved by semi-analytic approaches because of the requirement 
of the continuity in cosmological parameter space. It is thus desirable to understand the 
essential physics that dominates the dynamical evolution of intracluster gas. 

While intracluster gas is mainly driven by gravitational shocks and adiabatic compres- 
sion, two completely different mechanisms have been suggested thus far that may significantly 
affect the distribution of the intracluster gas: nongravitational preheating and radiative cool- 
ing. This is primarily motivated by the steepening of the X-ray luminosity-temperature rela- 
tion and the entropy excess of groups and clusters reported by a number of investigators over 
the past years (e.g. David et al. 1993; Wu, Xue & Fang 1999; Ponman, Cannon & Navarro 
1999). However, it has been realized that the prevailing preheating scenario is facing the 
so-called energy crisis - an unreasonably high efficiency of energy injection into the intraclus- 
ter medium from supernovae must be required (Wu, Fabian & Nulsen 2000). Moreover, a 
uniform preheating of the cosmic baryons to a temperature of ~ 10^ K would make the Lya 
forest to disappear at high redshifts. In contrast, radiative coohng is a natural process of the 
hot intracluster gas, in which no exotic physics is needed. Several recent studies have shown 
that radiative cooling alone can allow one to successfully reproduce not only the observed 
distributions of global X-ray luminosity and entropy but also the internal structures of hot 
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gas in groups and clusters (Bryan 2000; Pearce et al. 2000; Muanwong et al. 2001, 2002; 
Voit & Bryan 2001; Wu & Xue 2002a,b; Borgani et al. 2002; Voit et al. 2002). If confirmed, 
this would have profound implications for our understanding of the evolution of hot gas and 
the formation of stars in the most massive systems in the universe. 

So far, only can one employ hydrodynamical simulations to assess the effect of radiative 
cooling on the SZ counts and power spectrum (da Silva et al. 2000; White et al. 2002), while 
in the preheating scenario several (semi) analytic models have been applied for the prediction 
of SZ cluster surveys (Cavaliere & Menci 2001; Holder & Carlstrom 2001; Benson, Reichardt 
& Kamionkowski 2002). Now, a fully analytic treatment of radiative cooling effect on the 
prediction of SZ cluster counts and power spectrum may become possible if the observed 
central entropy floor of groups and clusters can be attributed to radiative cooling of hot 
intracluster/intragroup gas (Voit & Brayn 2001). In such a scenario, the minimum entropy 
of gas distribution in each cluster can be uniquely determined by the conservation of energy, 
which is equivalent to specifying the equation of state for the intracluster gas. Consequently, 
we will be able to derive the gas distribution of clusters at various redshifts as a result 
of cooling, in combination with the hydrostatic equilibrium hypothesis, provided that the 
underlying dark matter profile can be approximated by some kinds of analytic form, e.g., the 
universal density profile (Navarro, Frenk & White 1997; NFW). Together with the Press- 
Schechter (1994; PS) formahsm for the abundance of dark halos at different cosmic epoch, 
we can eventually compare the predicted SZ counts and power spectra with and without 
radiative cooling, and demonstrate the uncertainty in the determination of the cosmological 
parameters arising from radiative cooling. 



2. Gas distribution with and without cooling 

2.1. Dark matter distribution 

Intracluster gas with and without cooling is always assumed to be in hydrostatic equi- 
librium with the underlying gravitational potential dominated by dark matter component 

1 djn.kBT) _ G'MDM(r) 

J 9 ' \ ) 

IdmpHe dr r'^ 
where and T are the electron density and temperature, respectively, and = 0.585 is 
the mean molecular weight. We use the universal density profile suggested by numerical 
simulations (NFW) for pcM 

(r/rs)(l + r/rs)2 
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where 5c and are the characteristic density and length of the halo, respectively, and pcrit 
is the critical density of the universe at cosmic time t. We follow the prescription of Eke, 
Navarro & Steinmntz (2001) to fix the two free parameters, Sc and r^, in the NFW profile. 
To do this, the concentration parameter c = Tyu/^s of a dark halo identified at redshift z is 
related to the collapsing redshift Zcow through 

3 ^ Ac(^coll) ^m(^) / 1 + ^colA ^ /ON 

where Ac is the overdensity of dark halo within virial radius Tyir with respect to pcrit, for 
which we take Ac = IStt^ + 82[Qm(-2) — 1] — 39[Qm(-2) — 1]^ for a flat universe, and Qm(-2) is 
the cosmic density parameter. The collapsing redshift 2;coii is determined by 

^(^coll)(Teff(M,) = (4) 

(no- 
where D is the normahzed linear growth factor, = 25 for ACDM model, and aes{Ms) is 
the so-called modulated rms linear density at mass scale Mg 
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and Ms corresponds to the mass contained within r = 2.17rs where the circular velocity 
reaches the maximum. Finally, we specify a temperature T^i^ to the hot gas in cluster of 
mass M at redshift z in terms of cosmic virial theorem (Bryan & Norman 1998) 

^BT.i, = 1.39/T(/.^Aci?^)V3 keV [^-^^^J , (6) 

in which we will choose the normalization factor to be fr — 0.92, E'^ — J1m(1 + z)^ + + 
(1 — Qm — ^^a)(1 + and h is the Hubble constant in units of 100 km s~^ Mpc~^. 



2.2. Gas distribution without cooling 

In the absence of cooling, we adopt two models for the gas properties: 

Model I - Gas is assumed to follow the same distribution as the dark matter in clusters: 

/b /rr\ 

= Pdm (7) 

in which /b is the universal baryon fraction, and //e = 2/ {l-\-X) is the mean electron weight 
with X — 0.768 being the hydrogen mass fraction in the primordial abundance of hydrogen 
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and helium. We solve the equation of hydrostatic equilibrium to get the temperature profile 
under the boundary restriction T'^{r — > oo) 0: 

k.T ir) = k.T - ^1 + - j l^^^ .3(i + ,)3 d., (8) 



where ksT* — iiiG ^mpdcpcntf s- 

Model II - Gas is assumed to be isothermal and T^{r) — Tyir- In this case, the electron 

number density in terms of the equation of hydrostatic equilibrium reads (Makino, Sasaki & 
Suto 1998) 

n°(r)=neoe-Ml + -j , (9) 

where a — A-kG ^m^dcpciitfl/ k-eT . The normalization parameter can be fixed through 

AT:n^m^r^nl{r)dr = M^ir/b- (10) 



Jo 



2.3. Gas distribution with cooling 

The coohng time scale for a steady, highly subsonic flow is determined by the conser- 
vation of energy. Setting the energy loss rate due to bremsstrahlung emission to equal the 
change in the speciflc energy of gas yields 

— »'v(^)(^)'7^)-\ 

where and T° are the electron number density and temperature without cooling, respec- 
tively, and g is the total Guant factor. If the cooling time tc is chosen to be the age of 
clusters, or approximately the age of the universe, the above equation would allow us to set 
up a link between and T°. Note that this may lead to an overestimate of the cooling 
effect, and meanwhile tc becomes to be cosmological model dependent. Perhaps, a more rea- 
sonable approach is to deflne the age of a cluster as the cosmic time between its collapsing 
redshift Zcou and the redshift z when it is identifled (cf. Section 2.1). Once the gas density 
is specified, we will be able to work out the cooling radius and the amount of gas that cools 
out of the hot phase by the time t^, (e.g. Wu & Xuc 2002b). The fate of the cooled materials 
within cooling radius may be associated with star formation. Now, our task is to work out 
the new distribution of the intracluster gas with cooling under different assumptions of the 
equation of state. 
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Model III - Before and after cooling, the equation of state for gas is always assumed 
to be isothermal and T{r) — Tyij.. Consequently, the functional form of equation (9) also 
applies to the cooling case. However, the normalization after cooling is made through 

/ 47r/Xempr^ne(r)(ir = Myir/b - Mcooi (12) 

where the total cooled material Mcooi is given by 

Mcooi = / 47r^e"7.pn°(r)r^dr, (13) 
Jo 



= 47r/b5cPcritr 



1 I 1 , ^cool \ ^cool 

In IH 



(14) 



~l~ '"cool 

and Tcooi is the cooling radius which can be obtained by combining equations (9) and (11). 

Model IV - Gas traces dark matter (Model I) before coohng, and the equation of state for 
the new gas distribution after cooling is given by the entropy before cooling plus a constant 
entropy floor (cf. Holder & Carlstrom 2001; Voit & Bryan 2001) 

where the entropy floor is determined by the cooling time tc 

^B^vir 



= lOOkeVcm^f '^Y"(^f'(>^Y\ (16) 

V2.869 X 101V/ ^l-S/ V keV y ^ ^ 

This "entropy floor" increases with cosmic time, resulting in a deposition of cooled materials 
in the central regions of clusters (Voit & Bryan 2001). Solving the equation of hydrostatic 
equilibrium yields 

kBT{r) = ~G^m,S'/'{r) T ^^^S-'/'{r)dr, (17) 

^ Jr 

and ne = {kBT/Sf^. 

Model V - The same as Model IV except a varying metallicity Z = 0.3ZQ{t/to) is 
assumed instead of Z = 0.3^0 for the rest four models, where to denotes the present cosmic 
epoch. The parameters of the five models are summarized in Table 1. 
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3. Expectation for SZ cluster counts 

The total SZ flux observed at frequency u hy a cluster of mass M at redshift z is (e.g. 
Barbosa et al. 1996) 

S^(x, M, z) = (^) I k^T(r)n,(r) Anr'dr, (18) 

D^{z) \mecy J 

where 

. X ^ (faTcMB)^ x^e"" 
9^^) - 2^^-^^-^/.(x); (19) 

U{x) = xcoth|-4, (20) 

X — hpu/ksTcMB, TcMB — 2.728 K is the temperature of CMB (Fixsen et al. 1996), and Da 
is the angular diameter distance to the cluster. In the isothermal case, the integral in the 
right hand of equation (18) can be replaced by the total mass of the cluster if the universal 
baryon fraction is introduced: 

S. = ^(^](^) M. (21) 



This is the standard method adopted in the hterature for the theoretical prediction of SZ 
cluster counts (Model II). While the assumption of isothermality is more or less reasonable 
in terms of current X-ray observations, the dependence of the gas fraction on temperature 
claimed by many observations throws doubt upon the direct utilization of equation (21). 
Nevertheless, if the gas fraction fh is allowed to vary according to mass or temperature, or 
the total gas mass is replaced by the hot gas component alone, Mgas(T), we may modify the 
above equation to be 

S. = ^(^)(^)m,^(T,. (22) 



Table 1: The parameters and legend for the models of gas distribution. 



Model 


cooling 


rie 


T 


metallicity (Z©) 


line-style 


I 


no 


gas-traces-mass 


e.h.e.* 


0.3 


dot-dash 


II 


no 


e.h.e. 


isothermal 


0.3 


dashed 


III 


yes 


e.h.e. 


isothermal 


0.3 


dotted 


IV 


yes 


gas-traces-mass 


e.h.e. 


0.3 


sohd 


V 


yes 


gas-traces-mass 


e.h.e. 


0.3{t/to) 


dash dot dot dot 



*Obtained by solving the equation of hydrostatic equilibrium. 
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As shown in the above section, Mgas{T) can be analytically determined for any cluster within 
the framework of radiative cooling. This provides a simple approach to estimating the effect 
of radiative cooling on SZ cluster counts without knowing the density distribution of the hot, 
isothermal gas inside clusters, which will be applied to Model 111. For rest three models, I, 
IV and V, we will adopt the exact formula of S„, equation (18), in the theoretical prediction 
of SZ cluster counts. 

The expected number of SZ selected clusters with flux greater than S^, and in redshift 
interval {z, z + dz) is 

^^(> S.) r dn^^ 

dzdQ dzdfl J m^u^,s.) 
where the mass threshold M-^i^^z, S^) is given by the definition of equation (18) or (21) or 
(22), depending on what approximation we would use. We adopt the PS mass function to 
describe the distribution of clusters 

[2 p 6,{z) da f 6l{z)\ 
dn^ -a/-tt o ,.. cxp — ^ c^M, (24 

where p is the mean cosmic density, 6c is the linear over-density of perturbations that col- 
lapsed and virialized at redshift z, and a is the linear theory variance of the mass density 
fluctuation in sphere of mass M: 

oo 

a^{M) = 2^y" k''P{k)\W{kR)fdk, (25) 



and W{kR) — 3(sin a; — a; cos x)/x^ is the Fourier representation of the window function. The 
power spectrum, P{k) cx k"'T'^{k), is normalized by the rms fluctuation on an 8 h-'^ Mpc 
scale, (Tg, and we take the transfer function T{k) from an adiabatic CDM model given by 
Bardccn ct al. (1986) for the Harrison-Zel'dovich case n=l. Note that our results may be 
moderately changed if a modified mass function of dark halos is adopted (e.g. Jenkins et al. 
2001; Sheth & Tormen 2001). 

We work with a fiat cosmological model (ACDM) of Qbh"^ = 0.019, Qu = 0.3 and 
Qa = 0.7. We adopt cg = 0.90 and a Hubble constant oi h = 0.65. We choose the SZ fiux 
limit to be S„ = 15 mJy at frequency z/ = 30 GHz. We perform numerical calculations for 
the five models of gas distribution with and without cooling, and demonstrate the resulting 
differential SZ cluster counts in Figure 1. While there are some differences between the 
redshift distribution of clusters predicted by the gas-traces-mass assumption and the one by 
the isothermal assumption, the effect of radiative cooling leads to a moderate decrease in the 
expected number of SZ clusters although the peak locations with and without cooling remain 
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roughly the same. Alternatively, our prediction alters only slightly if the cosmic evolution 
of metallicity according to 0.3ZQ{t/to) is included. 

Recall that in the standard treatment, the expectation for SZ cluster counts is made by 
presumably taking the intracluster gas to be isothermal and without cooling. If this predic- 
tion is directly applied to the determination of cosmological parameters in future SZ surveys, 
large uncertainty may be introduced because of the ignorance of radiative cooling correction. 
Now we can evaluate the uncertainty in the determination of cosmological parameters by 
comparing the theoretical predictions with and without radiative cooling. To do this, we 
still work with a flat cosmological model of Qm + = Ij but allow to vary according to 
(jg until the predicted redshift distribution of SZ cluster counts by the standard model (II) 
matches those by the coohng models III-V. 

Figure 2 shows the resulting (Xg versus Qm, in which the contours indicate 68% joint 
confidence intervals on the two parameters, corresponding to Ax^ = 2.30. Uncertainties in 
flu and (78 due to radiative cooling correction can be easily demonstrated by noticing that 
the true input values in the cooling models are JIm = 0.30 and as — 0.90, while our best-fit 



results are (Qm,(J8) = (0.241°;^?, 0.891°:?^), (0.23l°:Ji, 0.88ino) and (0.231°;^?, 0.891°;°^) for 



the three coohng models III, IV, and V, respectively. Namely, without the correction of ra- 
diative cooling we might underestimate both and parameters for ACDM cosmological 
model, although within the uncertainties of our 'observations', the best-fit parameters are 
still consistent with the input values. 



The angular power spectrum of temperature fluctuation on the CMB sky due to the SZ 
effect of clusters can be separated into the Poisson term and clustering term (Cole 
& Kaiser 1988; Komatsu & Kitayama 1999): 



4. SZ Power spectrum 




(26) 



(27) 



and 




(28) 
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where z^ec ~ 1000 is the CMB photon decouphng redshift, Do is the comoving distance 
to cluster of mass M at z, yi{M,z) is the Fourier transform of the Compton y-parameter 
in the thermal SZ effect: AT/Tcmb = fv{.x)y{0), fu{x) reflects the spectral dependence, 
and b{M, z) is the so-called bias parameter, for which we use the analytic approximation of 
Mo & White (1996). In our numerical computation, the minimum cluster mass is taken to 
be Mmin = 1 X IO^^Mq, and our final results are unaffected by this choice if we focus on 
large-scale (>arcminute) fluctuations on the CMB sky. 

We illustrate in Figure 3 the SZ power spectra with and without cooling for the five 
models in Table 1, together with the primary CMB signal produced by CMBFAST. It appears 
that the power spectra predicted by the gas-traces-mass and isothermal assumptions become 
indistinguishable, and the total SZ power spectrum is dominated by the Poisson distribution 
of groups and clusters. Inclusion of radiative coohng leads to a significant decrease of SZ 
power spectrum, especially at small angular scales or high-£. This is simply because cooling 
removes more efficiently the hot gas at the central regions of clusters and in low-mass groups, 
and therefore suppresses the SZ signal at small scales. This result is in good agreement with 
the one found by hydrodynamical simulations (e.g. da Silva et al. 2001). 



5. Discussion and conclusions 

As a natural process, radiative cooling plays an important role in the formation of 
galaxies. In the central regions of clusters and groups, the typical cooling time scale of the 
hot gas is usually less than the age of the universe. It follows that a considerable amount of 
the intracluster/intragroup gas must have cooled out of the hot phase since the formation of 
clusters and groups unless the cooling efficiency can be suppressed by other non-gravitational 
heating processes such as energy feedback from supernovae and AGN activity. Since cooling 
removes the hot gas in groups and clusters, it meanwhile reduces the SZ signal, leading to 
a decease of the expected SZ cluster counts and power spectrum. This scenario has been 
confirmed by numerical simulations (e.g. da Silva et al. 2001) and our calculations based on 
a simple analytical model. If the excess entropy detected in the central cores of groups and 
clusters can be attributed to radiative cooling (Voit & Bryan 2001), some conclusions drawn 
from a phenomenological preheating model based on the observed entropy distribution can 
be equally applied to the cooling scenario. Indeed, it has been shown that the preheating 
model also results in a decrease of SZ cluster counts and power spectrum (e.g. Komatsu & 
Kitayama 1999; Holder & Carlstrom 2001; da Silva et al. 2001). Recall that current X-ray 
observations alone are still unable to distinguish between preheating and cooling models (e.g. 
Voit et al. 2002; Borgani et al. 2002). 
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lt is well known that radiative cooling is a runaway process. Consequently, the expected 
cooled materials within the framework of radiative cooling from both hydrodynamical sim- 
ulations and analytical approaches exceed the observed stellar mass fraction in the local 
universe (Balogh et al. 2001; Wu & Xue 2002b and references therein). Inclusion of star 
formation and energy feedback, which should in principle resolve the so-called cosmic cooling 
crisis, may significantly reduce the effect of radiative cooling on the X-ray properties and 
SZ signal of clusters. Indeed, recent numerical simulations by White et al. (2002) have 
shown that as a consequence of the combined effect of energy injection by star formation 
and radiative cooling, the SZ power spectrum remains roughly unchanged. 

Our simple analytic approach to estimating the effect of radiative cooling on the non- 
targeted SZ cluster survey and SZ power spectrum yields a result consistent with those found 
by hydrodynamical simulations (e.g. da Silva et al. 2001). This analytic method allows an 
effective evaluation of the significance of various factors in the gas cooling process. For 
example, we have shown that within the framework of radiative cooling, both SZ cluster 
counts and power spectrum are insensitive to temperature variation. In other words, the 
isothermality hypothesis is a good approximation in the theoretical study of cluster SZ effect. 
Also, one may choose to use the gas-traces-mass assumption, which provides a reasonable 
description of the gas distribution outside the X-ray core. Recall that unlike the X-ray 
emission which depends on n^, the SZ signal is proportional to rig. So, the gas outside the 
X-ray core makes an important contribution to SZ signal too. Finally, in order to eliminate 
the cooling crisis in the present model and include another natural process, i.e., heating 
by super novae, in the evolution of hot gas, it is worth exploring how energy injection into 
the intracluster/intrgroup gas from star formation is incorporated into cooling scenario in a 
simple analytic way. 

This work was supported by the National Science Foundation of China, and the Ministry 
of Science and Technology of China, under Grant No. NKBRSF G19990754. 
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z 



Fig. 1. — Expected rcdshift distribution of SZ cluster counts with Sy = 15mJy at frequency 
30 GHz. The results for five models with and without radiative cooling and under different 
assumptions of the intracluster gas in Table 1 are shown. Also illustrated in the inset is the 
minimum mass threshold against cluster redshift. 
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Fig. 2. — 68% joint confidence contours for Qm and obtained by minimizing tlie 
quantity: Tlie redsliift distribution of SZ cluster counts predicted by standard model (II) is 
required to match those by cooling models (see Table 1 for legend). A flat cosmological model 
of + = 1 is assumed. For the cooling models we have adopted (I^m? o"8) = (0.30, 0.90) 
(indicated by the filled circle). The cross, asterisk and plus symbols represent the best-fit 
results for Model III, IV and V, respectively. 
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Fig. 3. — SZ power spectra from intracluster gas with and without radiative coohng (see 
Table 1 for legend). Top and bottom curves show the Poisson and clustering contributions, 
respectively. The upper thick solid line is the primary CMB signal predicted by CMBFAST. 
Recent observational results from BIMA (open squares; Dawson et al. 2001) and CBI (filled 
circles; Mason et al. 2002) are also shown. 



